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ABSTRACT 

Magnetic dissipation is frequently invoked as a way of powering the observed emission 
of relativistic flows in Gamma Ray Bursts and Active Galactic Nuclei. Pulsar Wind 
Nebulae provide closer to home cosmic laboratories which can be used to test the 
hypothesis. To this end, we analyze the observational data on the spindown power 
of the Crab pulsar, energetics of the Crab nebula, and its magnetic field. We show 
that unless the magnetic inclination angle of the Crab pulsar is very close to 90 
degrees the overall magnetization of the striped wind after total dissipation of its 
stripes is significantly higher than that deduced in the Kennel-Coroniti model and 
recent axisymmetric simulations of Pulsar Wind Nebulae. On the other hand, higher 
wind magnetization is in conflict with the observed low magnetic field of the Crab 
nebula, unless it is subject to efficient dissipation inside the nebula as well. For the 
likely inclination angle of 45 degrees the data require magnetic dissipation on the 
timescale about 80 years, which is short compared to the life-time of the nebula but 
long compared to the time scale of Crab's gamma-ray flares. 

Key words: ISM: supernova remnants - MHD - magnetic fields - radiation mecha- 
nisms: non-thermal - relativity - pulsars: individual: Crab 



1 INTRODUCTION 

Magnetic fields are often invoked in models of the relativistic 
jet production by central engines of Active Galactic Nuclei 
(AGN) and Gamma Ray Bursts (GRB). In these theories 
the jets are Poynting-dominated at the origin with the mag- 
netization parameter a = B 2 /Anpc? 3> 1. This is different 
from the earlier essentially hydrodynamic, low a, models of 
relativistic jets in one important aspect. Even strong, high 
Mach number shocks, in high a plasma are weakly dissi- 
pative compared to their low a counterparts (e.g. Kenne 
& Coroniti||1984a| |Komissarov||2012l). Moreover, PIC sim 



ulations show that the acceleration of nonthermal particles 



may also be problematic at such shocks ( Sironi & Spitkovsky 



2009 



2011al. This suggests that either the Poynting flux is 
first converted into the bulk motion kinetic energy via ideal 
MHD mech anism (e.g. |Vlahakis fc Konigl 2004, Komissarov 



objects much "closer to home" which share similar proper- 
ties, the Pulsar Wind Nebulae (PWN). They are powered 
by relativistic winds from neutron stars and these winds are 
also expected to be Poynting-dominated at the base (see 
Arons|2012 and references therein). In particular, the Crab 
nebula, which is rightly considered as a testbed of relativis- 
tic astrophysics, has been studied with the level of detail, 
which may never be reached in studies of AGN and GRB 
jets. 

The early attempts to built a theoretical model of the 
Crab nebula using the ideal relativistic MHD approximation 
resulted in a paradoxical conclusion that the p ulsar wind has 
to have a ~ 10 -3 near its termination shock (Rees & Gunn 
|1974| |Kennel fc Coroniti||1984a| ). A slightly higher magne- 
tization, a ~ 10~ , was later suggested by axisymmetric 



et al.|[2009 Lyubarsky||2010 1 and then dissipated at shocks 
or the magnetic energy is converted directly into the energy 
of emitting particles via the magnetic dissipation which ac- 



companies magnetic reconnection events (e.g.|Drenkhahn & 


Spruit||2002| |Lyutikov & Blandford||2003| |Zhang 


& Yan 


2011 


Giannios|2011 


McKinney & Uzdensky|2012 


i. In fact, 



the magnetic dissipation can facilitate bulk acceleration of 
jets as well. 

While AGN and GRBs are very distant sources, which 
makes their observational studies more difficult, there exist 



numerical simulations ( Komissarov fc Lyubarsky||2003 Del 
|Zanna etal . 2004; B ogovalov et al.|2005[ ), although no proper 
study of this issue has been carried out. The key property of 
these analytical and numerical solutions is purely toroidal 
magnetic field. The strong hoop stress of such field creates 
excessive axial compression of the nebula in solutions with 
higher a and pushes the termination shock too close to the 
pulsar in the Kennel-Coroniti model, in conflict with the ob- 
servations. On the other hand, the ideal relativistic MHD 
acceleration of uncollimated wind-like flows is known to 
be very ineffective, leaving such flows Poynting-dominated 
on the astrophysically relevant scales (e.g. Lyubarsky|2011 
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Komissarov 2011 1. This striking conflict is known as the a- 
problem. 

Attempts have been made to see if a can be reduced 
via magnetic dissipation in the so-call striped zone of the 
pulsar wind, where the magnetic field changes it polarity 
on the length scale A p = cP, where P is the pulsar period 
( |Coroniti|1990l |Lyubarsky fc Kirk|2001| ). The dissipation is 
accompanied by the wind acceleration via conversion of the 
thermal energy into the bulk kinetic energy of the flow dur- 
ing its adiabatic expansion. Unfortunately, for the wind of 
the Crab pulsar the dissipation length scale significantly ex- 
ceeds the radius of the wind termination shock, thus making 
this mechanism inefficient ( Lyubarsky fc Kirk|2001 1. 



Lyubarsky ( 2003b I has demonstrated that the energy 



associated with the alternating component of magnetic field 
of the striped wind can be rapidly dissipated at the termi- 
nation shock itself, where the characteristic Larmor radius 
of shock-heated plasma exceeds the stripes wavelength. His 
solution of the shock equations, which accounts for the "eras- 
ing" of stripes, shows that the post-shock flow is the same 
as it would be if the dissipation had already been fully com- 



pleted in the wind. Sironi & Spitkovsky (2011b I have used 



3D PIC simulations to study the magnetic dissipation and 
particle acceleration at the termination shock of the striped 
wind numerically and concluded that efficient magnetic dis- 
sipation occurs even when the Larmor radius remains be- 
low the stripes wavelength, via rapid development of the 
tearing mode instability and magnetic reconnection in the 
post-shock flow. 

One way or another, this dissipations occurs only in the 
striped zone and only the alternating component of mag- 
netic field dissipates. Outside of the striped zone, around 
the poles, the pulsar wind a remains unaffected by this dis- 
sipation and hence very high. As we shall show, the overall 
magnetization of plasma injected into the nebula will also 
be rather high unless the striped zone spreads over almost 
the entire wind. 



Lyubarsky (2003a I argued that in the polar zone the 



wind a can be reduced via the flow acceleration related to 
the dissipation of fast magnetosonic waves emitted by the 
pulsar into the polar zone. However, it seems unlikely that 
the energy flux associated with these waves can dominate 
the wind energetics in the polar zone. At least, the 3D nu- 
merical simulations of pulsar winds (A. Spitkovsky, private 
communication) show that their contribution is rather small. 
Thus, we do not expect a of the polar zone to be below unity. 
An alternative solution to the a problem has been pro- 



posed by Begelman ( 1998 ) , who argued that the axial com- 



pression of the nebula can be reduced via the current-driven 
kink instability, resulting in more or less uniform total pres- 
sure distribution inside the nebula. This would make the 
overall structure and dynamics of the nebula similar to that 
in the models with particle-dominated pulsar wind. Recent 
computer simulations of the non-linear development of the 
kink instability of relativistic Z-pinch configurations support 
this conclusion ( Mizuno et al.|2009 20111. In this scenario, 
PWN are supplied with highly magnetized plasma, making 
magnetic dissipation a potentially important process in their 
evolution and emission. 

In this paper, we test whether the magnetic dissipation 
inside PWN is consistent with the observations of the Crab 
nebula and its pulsar. The main idea is very simple. First, 



the timing observations of the Crab pulsar allow us to esti- 
mate how much energy has being pumped into the nebula. 
Second, using the stripe wind model we can calculate how 
much of this energy is supplied in the magnetic form. Third, 
a simple dynamical model of the nebula expansion can be 
used to predict how much magnetic energy in retained by 
the nebula after adiabatic losses. Finally, the observations 
of the Crab nebula tell us how much magnetic energy is 
actually in there and whether the magnetic dissipation is 
actually required to make the ends meet. 



2 OVERALL ENERGETICS OF THE CRAB 
NEBULA 

In the simplest approximation, the spindown of pulsars is 
described by the equation Cl oc — f2™, where f2 is the pul- 
sar angular frequency and n is the so-called brake index. 
This form of the spindown law originates from the magneto- 
dipole vacuum radiation mechanism which gives n = 3. 
Force-free (or magnetodynamic) models of pulsar magneto- 
spheres yield the same dependence on Q ( Spitkovsky||2006 



Kalapotharakos & Contopoulos 20091. The solution to this 
equation is 



fi = fi 1 



The corresponding spindown luminosity is 

■ i + i 



-mn = L l + 



t 



(i) 



(2) 



where r is called the spindown time (Rees & Gunn 1974) 



From the timing observations of the Crab pulsar, n = 2.51 
and r ~ 703 yr ( Lyne et al.|1993 |. For the usually accepted 
moment of inertia of neutron stars I = 10 45 gcm 2 , the cur- 
rent spindown power L sp ~ 4.6 x 10 38 erg/s and the initial 



power L ~ 3.3 x 10 dy erg/s. 
energy of the Crab pulsar 



E = Lqt- 



l + ~ 



The total extracted rotational 



3.73 x 10 49 erg (3) 



or 67 percents of its initial rotational energy. The integrated 



luminosity of the Crab nebula L n ~ 1.3 x 10 erg/s (Hes- 



ter]|2008 1 is significantly below L sp suggesting that a large 
fraction of E is converted into the kinetic energy of the su- 
pernova shell and deposited as the internal energy of the 
PWN. In fact, both the optical and radio observations re- 
veal accelerated expansion of Crab's PWN and its network 



of thermal filaments ( Trimble 1968 Wyckoff fc Murray 1977 



Bietenholz et al.||1991[ |. The expansion velocity appears to 
have increased by 100-200 km/s during the lifetime of the 
nebula, which corresponds to the incr ease of the ki netic en- 
ergy of its filaments by E k ~ 10 49 erg ( |Hester||2008 1. 

Assuming that the magnetic field becomes randomized 
and behaves as gas with ultrarelativistic ratio of specific 
heats T = 4/3, the internal energy of PWN is E n = 3pV, 
where p is the PWN uniform total pressure and V is its 
volume. Ignoring the radiative cooling, the evolution of E n 
can be described by the equation 



E n — L s 



If v 

3 n V 



(4) 
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Figure 1. Left panel: Mean magnetization (cr a ) of the pulsar wind after dissipation of its stripes and the mean fraction (<5 Q ) of 
magnetic magnetic energy injected into PWN as functions of the pulsar magnetic inclination angle. Middle panel: The distribution of 
magnetic energy injected into PWN over the polar angle, f a (9) = (3/2)S a (9) sin 3 9, for the magnetic inclination angle a = 0° (solid line) 
15°, 30°, 45°, 60° (dashed lines). Right panel: The magnetization of the striped wind after dissipation of its stripes as a function of the 
polar angle for the magnetic inclination angle a = 0° (solid line), 15°, 30°, 45°, 60°, 75° (dashed lines). 



where the last term describes the adiabatic cooling. Since 
the Crab PWN shows only rather slow acceleration we may 
assume that its radius R oc t, which leads to 



IT' 



(5) 



Given the expression |2| for L sp , the initial condition 



E n (0) = 0, and assuming n < 3, 
this equation 



we find the solution to 



E n = 



Lqt 



3a + 2 x 



1 - 



(a — l)x 2 + ax + 1 
{x + 1) Q 



(6) 



where a — (n + l)/(n — 1) and x — t/r. For the parame- 
ters of the Crab pulsar this yields E n ~ 1.3 x 10 49 erg. The 
corresponding spindown energy converted into the kinetic 
energy of the supernova shell is then Ek — 2.4 x 10 49 erg, 
which agrees very well with the observational value given by 



Hester ( 2008 1 



3 THE MAGNETIC POWER OF STRIPED 
WIND 

If we ignore the magnetic dissipation (or amplification) in- 
side the nebula, as well as the radiative losses of particles, 
then the energy distribution between particles and magnetic 
fields in the nebula equals to that immediately downstream 
of the termination shock. In order to estimate the fraction 
of the wind energy remaining in the magnetic form we will 



employ the split-monopole model by |Bogovalov ( 1999 1 and 



the conclusion made in Lyubarsky (12003b I that the over- 



all effect of stripes dissipation at the termination shock is 
equivalent to their dissipation upstream of the shock. 

Let us denote as a the angle between the spin axis and 
the magnetic axis of the pulsar, the magnetic inclination 
angle, as 9 the angle between the rotation axis and the par- 
ticular streamline of the wind, and as <fi the phase of the 
stripe wave, <j> — corresponding to the middle of the stripe 
with positive (or negative) B^. Then the phases separating 
the positive and negative stripes are <f> a {9) and 2tt — <f>a (0) 



where 



cos0 a (#) = — cot(a) cot(6 l ). 



The conservation of the total magnetic flux corresponding to 
one wavelength allows us to find the magnitude of magnetic 
field after the dissipation of stripes as 



B = B 



|2<M#)A- 



7r/2 - a < 9 
6<n/2-a 



< tt/2 



where Bo is the magnitude of the magnetic field of the 
striped wind (In [Lyubarsky (2003b), B is called the mean 
magnetic field of the striped wind). In these calculations we 
assumed that after the dissipation of the alternating com- 
ponent of the magnetic field the relic current sheets col- 
lapse following their adiabatic cooling. The fraction of the 
wind power remaining in the form of Poynting flux along 
the stream line with the polar angle 9 is 



(2^(60/71-- 



1. 



7r/2 - a < 9 < tt/2 
6 < tt/2 - a 



Neglecting the small contribution of the bulk kinetic energy 
to the wind power before the stripes dissipation, the wind 
magnetization along the stream line after the dissipation is 



<r a {6) 



x4g) 
-X«(0) 



We define the mean magnetization of the wind, {u a ), as 
the ratio of its total Poynting flux to its total bulk kinetic 
energy flux. Since in the split monopole model the energy 
flux density varies with 9 as sin 2 9, 

(Xa> 



{(T a ) 



where 



(x4 



(8) 



tt/2 

(x«> = f Xc(0)i 



siii 9d9 . 



The mean magnetization is shown in the left panel of Fig- 
ure [1] One can see that unless the pulsar is almost an 
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Table 1. 



Q 


10° 


20° 


30° 


40° 


50° 


60° 


70° 


80° 


(So) 


0.82 


0.64 


0.48 


0.34 


0.23 


0.13 


0.061 


0.017 



orthogonal rotator its value is much higher compared to 
(er) ~ 10 -3 of the Kennel-Coroniti model and the va lues uti- 
lized in the 2D numerical simulations, (a) ~ 10 -2 (Komis- 



sarov fc Lyubarsky||2003| |Del Zanna et aT7||2004| |Bogovalov 



et al.||2005 Camus et al. |2009 1. Unfortunately, a is poorly 



constrained from observations. Using as a guide the value 
obtained from fitting the spectrum and pulse profile of the 
high energy emission of the Crab pulsar, a ~ 45° ( |Harding 
|et al.|[2008[ l, we obtain (a) ~ 0.26. Thus, the dissipation 
of magnetic stripes is apparently unable to resolve the a- 
problem completely. The left panel of Figure [I] shows the 
distribution of a over the polar angle, where its value out- 
side of the striped zone is artificially limited by the rather 
arbitrary value of ~ 100. In reality, its value there could be 
determined by the dissipation of fast magnetosonic waves 
emitted by the pulsar (Lyubarsky 2003a I. However, the effi- 
ciency of this emission in 3D numerical simulations of dipo- 
lar pulsar magnetospheres seems to be rather low and thus 
one would indeed expect a rather high magnetization in the 
polar region. 

Next we consider the compression of the plasma of such 
a wind at the termination shock. The magnetic flux con- 
servation ensures that at the shock Bv n — const, where 
v n is the normal component of velocity. This implies that 
the Poynting flux increases by the shock compression factor 
T) — «n,i/fn,2- In the case of strong ultrarelativistic shock, 



■nix) = 6 (l + X + Vl + Mx + X 2 ) 



(9) 



This result holds not only for a perpendicular shock but also 
for an oblique shock (see Eq.A14 in |Komissarov fc Lyutikovj 
(2011 1). Thus, the fraction of energy injected into PWN in 



the magnetic form along the given streamline is 



(10) 



In the split monopole model the overall fraction of the wind 
power injected into PWN in the magnetic form is given by 
the integral 



tt/2 



<<U = 



5 a (0) sin 3 8d6. 



(11) 



The function (S a ) is shown in the left panel of Figure [T] and 
m Table [1] One can see that, unless the magnetic inclina- 
tion is close to 90° the fraction of magnetic energy is quite 
substantial. For the guide value of a ~ 45° (Harding et al. 
2008[), we obtain {6} ~ 0.28. Thus, almost one third of the 



energy supplied into the Crab nebula could be in the mag- 
netic form. The middle panel of Figure 1 shows how this 
flux is distributed over the polar angle for different mag- 
netic inclinations. For a < 50° it peaks at the boundary of 
the striped zone, but for a > 50° the maximum is inside the 
striped zone. 



4 MAGNETIC DISSIPATION INSIDE THE 
NEBULA 

The observed synchrotron and inverse-Compton emission of 
the nebula is well fitted by the "one-zone" model with mag- 
netic field of strength B ~ 125 ( |Meyer et aLl|2010[ ). Al- 
though the magnetic field in the nebula is unlikely to be 
uniform, this estimate has to be more reliable compared to 
the usual equipartition one, which requires an additional as- 
sumption of parity between the energies of magnetic field 
and synchrotron electrons. The observed shape of the neb- 
ula can be described as a prolate spheroid with the major 
and minor axes a = 4.4 pc and b = 2.9 pc (Hester 20081, 
which gives the volume V = (Tv/6)ab 2 ~ 5.7 x 10 cm . 
The corresponding total magnetic energy of the nebula is 
E m = 3.5 x 10 47 erg, which is significantly below the value 
of E n we estimated in Sec [2] 

Assuming parity between the energy of synchrotron 
electrons E e and the magnetic energy E m , |Hillas et al.| 



( 1998 1 used the observed synchrotron luminosity of the neb- 
ula to derived its equipartition magnetic field B eq = 330 pG. 



The lower value of B given by Meyer et al. (20101 suggests 
significant deviation from the energy equipartition. From the 
theory of synchrotron emission it follows that the energy of 
electrons E e oc B -3//2 whereas E m oc B 2 . Thus 



E e = E m (B/B e 



-7/2 



In the case of the Crab nebula this yields E e ~ 30i5 m ~ 
1.0 x 10 49 erg, which is remarkably close to our value of E n . 

This result suggests two possible explanations. First, 
the energy may be supplied into the nebula mainly in the 
form of relativistic particles. The analysis presented in the 
previous section shows that this would require the Crab pul- 
sar to be almost an orthogonal rotator, in fact we would 
need a ~ 76°. If however the magnetic inclination angle is 
indeed close to a = 45°, obtained in Harding et al. ( 2008 1 via 



modeling of the pulsed emission, then efficient dissipation of 
magnetic field inside the nebula, accompanied by particle 
acceleration, is required to explain the data. 

Assuming that a fraction (8) of L 3p is supplied into the 
nebula in the magnetic form, one can find the characteristic 
timescale of this dissipation via balancing the supply and 
dissipation rates as 



En 



Trad — 



{S)L a p 



(A 

0.3 



yr. 



(12) 



This is much smaller compared to the dynamical timescale 
Tdn — 950 yr, which justifies the omission of adiabatic en- 
ergy losses in this estimate. Moreover, r m d exceeds the light 
crossing time of the nebula, ti c — 12 yr, only by a factor of 
~ 7. This implies multiple reconnection regions inside the 
nebula. Indeed, the speed of magnetic energy supply into 
the reconnection zone is likely to be limited from above by 
~ 0.1 of the Alfv'en speed ([Lyubarsky 2005), which in the 
relativistic MHD is 

1/2 



where a = B 2 /4-kw, where w = pc 2 + Tp/ (F — 1) is the rel- 
ativistic enthalpy and p is the gas pressure. In magnetically 
dominated plasma a 3> 1 and c a is close to the speed of light, 
whereas in particle-dominated plasma with a <C 1, it can be 



© 0000 RAS, MNRAS 000, 000-000 



Magnetic dissipation in the Crab nebula 5 



significantly lower. The mean a of the nebula can be esti- 
mated as <a>~ 2E m /E e ~ 0.07 leading to <c a >~ 0.25c. 
If the reconnection flow was regular on the scale of the whole 
nebula, like for example in Lyutikov (2010|, the correspond- 
ing dissipation timescale would be ~ 300 yr, which is sig- 
nificantly higher than that given by Eqjl2] Thus smaller 
magnetic domains must be involved. 

In fact, the high-resolution photometric observations of 
the Crab nebula do show significant structure in the form 
of arcs and filaments (e.g. |Bietenholz et aL|2004 k More- 
over, the polarimetric observations by |Bietenholz fc Kron-| 
berg ( 1991 1 reveal the polarization which in the central part 



5 DISCUSSION 

The 2D RMHD numerical simulations of the Crab nebula 



of the nebula is about 5 times below that for the synchrotron 
emission in uniform magnetic field. They explain this result 
by the presence along the line of sight of several cells with 
randomly oriented uniform magnetic field. 

There is a number of different factors that can leads 
to such substructure. For example, 2D axisymmetric simu- 
lations of the nebula by |Camus et al.| ( |2009[ > reveal strong 
variability of the termination shock and generation of inho- 
mogeneities with the characteristic length scale similar to 
the shock radius, which is only about several light months 
long. It can also be produced via dynamic interaction be- 
tween the relativistic flow and the massive filaments of line- 
emitting thermal plasma apparently threading the PWN. 
This is supported by the observed spatial correspondence 
between these filaments and the nonthermal radio filaments 
( jVelusamy et al.|1992| >. 

Finally, it could result from the development of the cur- 



rent driven kink instability near the polar axis (Begelman 



1998). As we have seen, the structure of striped wind im- 



plies supply of low a plasma with slower reconnection speed 
into the equatorial region and high a plasma with higher 
reconnection speed into the polar region of the nebula. In 
this polar region, the Chandra and HST observations reveal 
the jet-like flow, which is likely to be a plume of highly mag- 
netized plasma compressed and squeezed away by the mag- 



netic hoop stress ( Lyubarsky 2002 Komissarov & Lyubarsky 



2003||Del Zanna et al.|2004||Bogovalov et al.|2005[ ). This jet 
continues up to ~ 1 ly and then merges with the background. 
The pronounced kink near the end of the jet suggests its de- 
struction by the instability. 

Thus, it seems quite likely that the Crab nebula does 
indeed have magnetic field with space variation on the scale 
< 1 ly, required to speed up the magnetic dissipation. The 
characteristic timescale of magnetic reconnection in such 
compact high a regions ~ 10 yr, which is short compared 
to the magnetic dissipation timescale of Eq |12| This is ex- 
pected as the reconnection and the dissipation timescales 
can be the same only in the case of oppositely directed re- 
connecting magnetic fields, whereas in the presence of strong 
guide field only a fraction of the incoming Poynting flux is 
dissipated during the reconnection (Lyubarsky 2005). Thus, 
one would expect numerous reconnection events to occur be- 
fore the magnetic energy gets dissipated, though at present 
it is difficult to give a quantitative estimate and further in- 
vestigation is required. 



| Komissarov fc Lyubarsky|2003| [Del Zanna et al.|2004| [Bo 



govalov et al.|2005 Camus et al.|2009 \ have been very suc- 
cessful in reproducing many key properties of the nebula, 
such as its jet-torus, the brightness asymmetry, wisps, and 
even the bright "inner knot" (Hester et al.||1995). In agree- 



ment with the observations, the proper motion of the sim- 
ulated jet and wisps is relatively low, v = 0.2 — 0.7c, as 
expected downstream of an almost purely hydrodynamical 
shock wave. This success leaves little doubt that the models 
capture the physics of the nebula quite well. 

However, as we can see now, the overall low wind mag- 
netization utilized in these models, <er>~ 10~ 2 , is in con- 
flict with what we would expect in the striped wind model 
without imposing very large magnetic inclination angle of 
the pulsar. Retrospectively, no attempts have been made 
to study models with much higher a and hence one can- 
not claim yet that the numerical simulations rule them out. 
The choice of a has been influenced by the very low value 
required in the Kennel-Coroniti model in order to have a 
termination shock in their ID solution. However, the flow 
dynamics of the 2D numerical solutions is completely differ- 
ent, as it involves large scale circulation and mixing, and we 
need to address this issue afresh. In fact, < a >~ 10~ 2 is 
already one order of magnitude higher compared to that of 
the Kennel-Coroniti modeQ Moreover, 3D solutions are ex- 
pected to be even more contrasting with the over-simplified 
ID solution. 



In Komissarov & Lyubarsky (20041, it was found that 



the size of the termination shock was somewhat decreasing 
with <<r>. As the shock size is determined by the balance 
between the wind ram pressure and the total pressure in the 
nebula, this tendency can be explained by the increasing ax- 
ial compression of the nebula by the magnetic hoop stress. 
However, this compression is certainly excessive in 2D mod- 
els, being enforced by the condition of axial symmetry which 
does not allow development of the kink instability. Thus, the 
ultimate answer to the question whether <a> 2> 10~ 2 are 
allowed by the RMHD model will only be found in future 
3D simulations. 

If at high latitudes the wind of the Crab pulsar is indeed 
free from stripes and has high a then downstream of the ter- 
mination shock one would expect a very fast flow with the 
Lorentz factor, 7 ~ a 1 ^ 2 in the case of perpendicular shock 



(Kennel & Coroniti 1984a I and even higher in the case of 
oblique shock ( Komissarov fe Lyutikov]|2011 1. Downstream 
of a perpendicular shock the flow is subsonic (or sub-fast- 
magnetosonic to be more precise), and can smoothly de- 
celerate down to 7 ~ 1 inside the nebula. Downstream of 
an oblique shock it may remain supersonic and a secondary 
shock will have to appear somewhere on its way. So far the 
observations of the Crab nebula show no evidence of such 
a secondary shock or such a fast flow. This may well be 
related to the low dissipation efficiency of shocks in highly 
magnetized plasma (e.g. Kennel fe Coroniti|| 1984a Komis 
sarov|20l"2 \ as well as the inability of such shocks to acceler- 



1 Komissarov fc Lyubarsky] | |2004| found that for < a >~ 10 -3 of 
the Kennel-Coroniti model the solution did not exhibit the polar 
jet. 
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ate non-thermal particles ( Sironi & Spitkovsky 2009 2011a I 



Further investigation is required to clarify this issue. 

The magnetic dissipation itself may produce visible and 
variable features in the images of the nebula, like arcs and 
filaments different from the line-emitting thermal filaments 



of the nebula. |Velusamy et al. ( 1992 1 state that VLA im- 
ages show such distinct features in the outer parts of the 
nebula. Given the very large synchrotron life-time of radio 
electrons compared to the age of the nebula one would not 
expect any significant difference in the radio spectral index 
of these arches and loops compared to the diffuse emission 
of the main body of the nebula, in general agreement with 
the observations (Bietenholz et al. 19971. This shows that 



radio observations are not really suitable for differentiating 
acceleration cites from cites with enhanced emissivity due 
to other reasons, e.g. stronger magnetic field. In the X-ray 
band, the life-time is much shorter and the spectral harden- 
ing could give away the ongoing particle acceleration. The 
X-ray observations of arcs and filaments in the outskirts of 
the Crab nebula do not show such hardening ( |Seward|2006[ ) , 
possibly ruling them out as the potential acceleration cites. 

Most of the radio filaments, particular in the inner part 
of the nebula, seem to coincide with the line-emitting fila- 
ments, though the author has been unable to find any recent 
work with quantitative analysis of this issue. Interestingly, 
these filaments do not seem to show up in the maps of polar- 
ized optical emission, which one would expect to see if the 
enhanced radio emissivity was simply due to the increase 
of magnetic field near the thermal filaments caused by the 
magnetic draping. This result may be important for the is- 
sue of the origin of radio electrons of the nebula. It is well 
known that the observed number of radio electrons is diffi- 
cult to explain in the current models of pair production in 
pulsar magnetospheres ( Arons|2012| ). However, this number 
is tiny compared to what is available in the line-emitting 
filaments. 

The magnetic reconnection/dissipation paradigm has 
received a boost by the recent observations of strong flares 



of gamma-ray emission from the Crab nebula ( Tavani et al 
2011 Abdo et al.|201 1 1 . However, the dissipation time scale 
given by Eq |12| is at least three orders of magnitude longer 
than the typical flare duration. It is possible that the tear- 
ing instability produces much smaller structures inside the 
large scale current sheets, however in this case one would 
expect a whole spectrum of time scales to be present. More- 
over, the current reconnection models of these flares involve 



strong magnetic fields, of order 1000 /iG (Uzdensky et al 



2011 Cerutti et al. 20121 and/or large bulk Lorentz fac- 



tors r > few ( Komissarov fc Lyutikov|2011 Clausen-Brown 



& Lyutikov 2012 1. Such conditions are not typical for the 



Crab nebula but, as we have pointed out, they are expected 
in the polar region near the termination shock, where the 
freshly supplied plasma can have very high magnetization 
and stream with ultra-relativistic speeds. Large Lorentz fac- 
tors could also be produced during fast reconnection events 
inside high-cr plasma, which again points out towards the in- 
ner polar region of the Crab nebula, where the observations 
reveal the Crab jet. It may well be that this is the most active 
region where most of the magnetic dissipation/reconnection 
takes place. 

The inner X-ray ring ( Hester et aL|2002 1 is another sus- 
pect. This is the most compact circular feature in the Crab 



nebula almost centered on the pulsar and often identified 
with the termination shock itself. Although the termination 
shock is not a ring-like feature but a closed surface surround- 
ing the pulsar, it may appear as a ring or a belt, particularly 
if the particle acceleration is confined to its equatorial sec- 
tion corresponding to the striped zone of the pulsar wind. 
The inner ring is a rather mysterious feature. It is not really 
seen in the optical images of the inner Crab nebula obtained 
with comparable angular resolution, apart from few fine de- 
tails ( jHester et al.|2002| . It is rather symmetric, defying the 
Doppler beaming effect which explains rather well the over- 
all asymmetry of the inner nebula within its MHD models 
(e.g. Komissarov & Lyubarsky 20031. The apparent symme- 
try of the ring is mainly due to its highly variably knots, 
rather than smooth underlying emission. The observed flar- 
ing of the knots makes them attractive candidates for the 
cites of magnetic dissipation, yet the fact that the ring is 
more or less centered on the pulsar suggests that it is lo- 
cated near the equatorial plane. 



6 SUMMARY 

(i) We have calculated the power of high-a striped pulsar 
wind which remains as the Poynting flux after total dissipa- 
tion of its stripes. The results show that the pulsar has to be 
an almost exact orthogonal rotator for the mean wind a to 
reduce down to the very low values suggested by the Kennel- 
Coroniti model and by the current axisymmetric numerical 
models of the Crab nebula. For the more realistic magnetic 
inclination angle a = 45°, about 30 percent of the wind 
power is retained in the form of the Poynting flux. While 
low magnetization is achieved in the equatorial plane, in the 
polar zone the magnetization remains very high. 

(ii) Given the relatively long spindown time of the Crab 
pulsar, and ignoring the radiative losses, we find that out of 
E ~ 3.7 x 10 49 erg, which have been supplied by the pulsar 
wind into the nebula, E n ~ 1.3 x 10 49 erg should still remain 
as its internal energy, sheared between magnetic field and 
relativistic particles. 

(iii) The observations of the synchrotron and inverse- 
Compton emission of the Crab nebula indicate that most 
of E n is stored in relativistic electrons and positrons and 
only Em ~ 3.5 x 10 47 erg in the magnetic field. This may 
be simply down to the fact that from the start the energy 
is injected into the nebula mostly in the form of relativis- 
tic particles. In the striped wind model, this implies that 
the Crab pulsar is almost an exact orthogonal rotator. Al- 
ternatively, most of the injected magnetic energy may have 
been dissipated and transferred to the particles via magnetic 
reconnection events. 

(iv) Using the magnetic inclination angle of the Crab pul- 
sar derived from modeling of its high energy pulsed emis- 
sion, a = 45°, we estimate the characteristic timescale of 
magnetic dissipation in the Crab nebula to be r m d ~ 80 yr. 
This relatively short timescale implies complex structure in 
the magnetic field distribution inside the nebula, which is 
supported by the radio and optical polarization data. 

(v) A combination of 3D magnetic instabilities and the re- 
sulting efficient magnetic dissipation may hold the key to the 
solution of the a-problem of the Crab nebula. Firstly, the 3D 
instabilities destroy the highly organized magnetic structure 
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of the oversimplified one-dimensional analytical solution of 
Kennel-Coroniti, and current axisymmetric numerical solu- 
tions, and thus reduce the axial compression of the nebula. 
This alone may be sufficient to allow much higher magneti- 
zation of the plasma injected into the nebula by the pulsar 
wind. Secondly, they create conditions which are favorable 
for magnetic dissipation, which reduce the dynamical effects 
specific to the magnetic field even further. 3D numerical sim- 
ulations are required to test this hypothesis. 

(vi) Since the scale of deduced magnetic dissipation inside 
the Crab nebula depends on the magnetic inclination angle 
of the Crab pulsar, accurate observational measurements of 
this angle are required in order to test its significance. The 
proximity of the nebula, makes possible search for in situ 
particle acceleration cites, but so far there have been no re- 
ports of clear signatures of ongoing magnetic reconnection 
in high resolution observations of the nebula. The recently 
discovered gamma-ray flares may be the first strong indica- 
tion of such activity. Further theoretical studies are needed 
to check if their short duration, only few days, is consistent 
with the much longer time scale of magnetic dissipation de- 
duce here. Such studies of the Crab nebula and other PWN 
will have important implications for the theory of AGN and 
GRB, which are much further out and cannot be studied to 
the same degree of detail. 
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